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Abstract 

In this paper, we propose a new approach to calculate multi-soliton solutions of Camassa- 
Holm (CH) equation and modified Camassa-Holm (MCH) equation with aid of Darboux 
transformation (DT). The new approach simplifies the approach presented in Proc. R. Soc. 
Land. A 460 2617-2627 (2004). We first map the CH and MCH equation to a negative order 
KdV (NKdV) equation by a reciprocal transformation. Then we proceed to apply the DT to 
solve the NKdV equation in the usual way. Finally we invert the reciprocal transformation 
to recover the solutions of the CH equation and MCH equation. 
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1 Introduction 

In recent years, the Camassa-Holm (CH) equation [T] 

qt + 2 qux + qxu = 0 , q = m + = u - u^x + , (1) 

{k is an arbitrary constant) derived by Camassa and Holm [T] as a shallow water wave model, 
has attracted much attention and various studies. This equation first appeared in the paper 
of Fuchssteiner and Fokas [2] on hereditary symmetries, but it came to be remarkable in the 
work of Camassa and Holm [T] where the peakon was described. As an integrable system, the 
CH equation shares most of the important properties of an integrable system of KdV type. For 
example, it admits Lax representation [T], bi-Hamiltonian formula mi, and can be solved by 
the inverse scattering transformation ii- The most remarkable feature of the CH equation 
is that it admits peaked soliton (peakon) solutions in the case A: = 0 m- ^ peakon is a weak 
solution in some Sobolev space with corner at its crest. The stability of the peakon were proved 
by Constantin and Strauss in the references ii- 

The interesting features of the CH equation have excited many endeavors to seek smooth 
soliton solutions of the CH equation. By a loop group approach, Schiff [S] found a single-soliton 
solution in a parametric form and an explicit (but incomplete) expression for the two-soliton 
solution. In the reference [5], Constantin presented a complete description of the scattering 
and inverse scattering problems associated with the CH equation. The key observation made 
by Constantin is that the isospectral problem of CH equation was mapped onto a Schrodinger 
spectral problem by introducing the Liouville transformation. Latter, Johnson [9] implemented 
the inverse scattering method for CH equation developed by Constantin to derive the soliton 
solutions of CH equation. Johnson succeeded in finding single-, two- and three-soliton solutions, 
but encountered difficulties in seeking the multi-soliton of the CH equation, since his procedure 
relies on a good guesswork and an extensive use of mathematical packages. To overcome the 
difficulties, Parker presented a different route to derive multi-soliton solutions of 

the CH equation based on Hirota’s bilinear transformation method. Parker first showed that 
the CH equation can be mapped onto a version of the AKNS-SWW equation HU by using the 
reciprocal transformation proposed in mm- Then Parker solved the AKNS-SWW equation 
by Hirota’s bilinear transformation method. Finally by the inverse reciprocal transformation, 
Parker arrived at the desired multi-soliton solutions of the CH equation. 

In this paper, inspired by the Parker’s work mmm, we propose an approach to compute 
the multi-soliton solution of CH equation with aid of Darboux transformation (DT). We remark 
that in m, Li and Zhang have refined Johnson’s method [9] for calculating explicit soliton 
solutions of the Camassa-Holm equation via DT. They found an analytical formula for the 
potential q and presented an algorithm for solving for the potential u. The Li-Zhang’s approach 
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is interesting and instructive to us, but the approach did not give an explicit formula for the 
multi-soliton solution for the potential u, since the algorithm requires one to solve second-order 
ordinary differential equations (ODEs) to compute the solution u. In contrast with Li-Zhang’s 
approach, our approach provides a more direct and simple way (avoid having to solve the ODEs 
encountered by Li and Zhang) to calculate the multi-soliton solution for the potential u of CH 
equation. We first make a connection between the CH equation and a negative order KdV 
(NKdV) equation by the reciprocal transformation proposed in [HI [16]. Then we proceed to 
apply the DT to solve the NKdV equation in the usual way. Based on the results of the DT, 
we present a procedure to invert the reciprocal transformation to recover the solution of the 
potential u. 

In this paper, we also show our approach can be applied to derive the multi-soliton solution 
of the modified CH (MCH) equation with cubic nonlinearity [HI [151 IISI EZl IB] 

mt -b [m{u^ - ul)] ^ = m = u-Uxx- (2) 

Similar with the case of the CH equation, we first map the MCH equation onto a NKdV equation 
by a reciprocal transformation. Then by using our procedure, we arrive at the multi-soliton 
solution of MCH equation. 

The whole paper is organized as follows. In Section 2, we shortly summarize the reciprocal 
transformation of CH equation, and describe our approach for computing the multi-soliton 
solution of CH equation. Then in Section 3 and Section 4, we apply our approach to calculate 
the multi-soliton solutions of CH equation and MCH equation, respectively. Some conclusions 
and discussions are addressed in Section 5. 

2 A reciprocal transformation for CH equation 

In this section, we first review some basic results about the reciprocal transformation of CH 
equation presented in [SIEKISI. These results enable us to relate CH equation to a NKdV 
equation. Then we propose our approach to calculate the multi-soliton solution of CH equation 
with aid of DT at the end of this section. 

It was shown in [T] that the CH equation possesses the following Lax pair with the spatial 
part 

il^xx = (^ + (3) 

and the temporal part 


( 4 ) 
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where A is a spectral parameter. 

The Liouville transformation [5] 

r = yjq = \/m + k‘^, (5) 

converts the CH equation into conservation form 

rt + {ur)x = 0. (6) 

This permits us to define a reciprocal transformation {x,t) (?/, t) by the relation 

dy = rdx — urdt, dr = dt, (7) 

and we have 

d d d d d 

di ^ ^ ^ ^ 

This reciprocal transformation was originally suggested in m and was later taken up in mum 

and mm- 

Under the reciprocal transformation ([7D, the CH equation ([T]) is transformed to a coupled 
system 

Xr + r'^Uy = 0 , 

u = r‘^ — rdy^ In r — , 

which was called as associated CH (ACH) equation [U |20] . 

Denote (f) = By using the reciprocal transformation ([7]), the spatial part of the spectral 
problem ([3]) becomes a Schrodinger spectral problem 

(f>yy = X(f)+U(f>, ( 11 ) 


(9) 

( 10 ) 


with the potential 


jj _ ^ 

2 r 


1 r 1 
_^ 

4 ; 7 *^ 


( 12 ) 


Now Li-Zhang’s procedure m for solving the CH equation using the DT is apparent. They first 
construct the potential U by applying the DT to the Schrodinger spectral problem (llip . Then 
they solve (1121) for r with the known U and solve u — Uxx + for u with the known r. 

By integrating ^ = x, they build up the relation between x and y with the known r. Finally, 
they introduce a variable t in the solution u of the CH equation. The Li-Zhang’s procedure is 
interesting and very instructive to us, but the difficulty encountered in executing the procedure 
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is that one have to solve the ODE (jl2p and the ODE u — Uxx + k‘^ = to complete the solution 
u. Inspired by the Parker’s work in the following we will accomplish the solution u of CH 


equation using DT by a quite different procedure, which enables us to avoid solving the ODEs 
encountered in Li-Zhang’s procedure. 

We notice that the reciprocal transformation ([7|) casts the temporal part of the spectral 
problem ([1]) into 



(13) 


The compatibility condition of (fTTI) and (fT^ gives rise to 



(14) 


We point out that (fTlt) is nothing but a negative order KdV (NKdV) equation [T91I201 [211 l22ll^ . 

Armed with the above results, we are now able to present a simple approach to derive the 
multi-soliton solution of CH equation via DT. The procedure is constituted of the following 
three steps: 

• apply Darboux transformation to NKdV equation (I14p to derive the solutions U and r; 

• substitute the known r into the formula (jlOp to calculate the potential tt; 

• build up a relation between the variables x, t and y, r. 

3 Application to solving CH equation 

Now we execute the procedure we presented above to calculate the multi-soliton solutions of CH 
equation. 

Step 1: Darboux transformation of the NKdV equation (llTh 

Let (po = <po{y,T) be a solution of Lax pair (ITTp and (IT^ fixed at the point A = Aq, and use it 
to define the gauge transformation 


(j) = {dy- a)(t), 


(15) 


where a = dy In cpQ. Let us set 

U = U — ay, f = r — Ur- (16) 


(16) 


By a direct calculation, we may check that (for the details see [23]) the Lax pair (flTT) and (ITSP 
is covariant with respect to the action of the Darboux transformation 


( 17 ) 
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Thus we arrive at the following proposition. 

Proposition 1 A solution U, r of the NKdV equation is mapped onto its new solution U, 
f under the Darhoux transformation m- 


Now we apply the above Darboux transformation to derive the multi-soliton solutions for the 
NKdV equation ([HI). Let /i, /2, •'' > /a^ be solutions of ifTT]) and (fT3]l fixed at A = Ai, A 2 , • • • , Xn- 
We start with a seed solution r = k and U = ^ and introduce the parameters pj, Cj by 
Ipj = Xj + Cj = • Then the functions /i, /2, • • • Jn can be chosen as 


/. 


coshf^j, j = 21 — 1, 
sinhf^j, j = 21, 


(18) 


where = ^pj{y+CjT). We define the Wronskian determinant Wn of N functions /i, /2, •'' > /at 
by 

WN = detA, Aij = , hi = ,V. (19) 

By repeating the Darboux transformation V-times (see the reference [H] for example), the 
V-soliton solutions of the NKdV equation (fHj) can be expressed as 

U = ( 20 ) 

r = k-2-^lnW(h,h,---jN)Ak-2^\nWn. (21) 

The two fundamental solutions of equation m with A = 0 and U given by (I2UI) are 


{ j, _ A WN 

W(h,f2,-,fN)^ - Wn^ 

X _ W{fi,f2 ,--- A 

W{h,f2,-,fN) Wn • 

Their asymptotic properties are 



N 

1 

2k 



N 

1 PA 

2k 2 ^ 


(/>! ~ n ( 

Pj 

2 . 

)e^y, (p2 • 

- n(- 

)e~^y, 


j-1 




i-1 




1 

~ 2k 

N 

■ 



N 

:-2i- 

4>l,y 

n( 

- ^)e^y, 

02,2/ ~ 

-^ni 







j-i 



( 22 ) 


(23) 


Step 2: the solution u{y, r) of CH equation 

The formula (j2ip provides a multi-soliton solution of the ACH equation (|9|). Substituting (j2ip 
into (fTOl) yields 

u = ^[{lnWN)yr? - m^WN)yr - [k - 2(ln Wjv)yr][ln(A: - 2(lnIT^)^,)]^,, (24) 

which provides an analytic solution of the CH equation expressed in the variables y and r. 
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Step 3: the relation between x, t and y, r 


To complete the solution u{x,t) of the CH equation, we have only to obtain the coordinate 
transformation between x, t and y, r. There are a number of ways to accomplish this final step 
(for example, see [U [9l UHl El])- We hnish this step based on a result of the DT presented in 
Step 1. When A = 0, equations (l3|) and m become 




(25) 

(26) 


4^yy — Ucj). 

For equation ([25]), we may choose a solution ip to he t/j = . For equation (f26ll . the solution cp is 

a linear combination of the two fundamental solutions (pi and (p 2 given by (I22p . The asymptotic 
behaviours (f23t) imply that the (p corresponding to ip = e "2 must be dicpi^ where di is a constant. 
Thus the relation (p = r^ip yields 

dicpi 


62 = 




which gives rise to the following relation 

dlcPl 


X = In ■ 


= In ■ 




+ cr, 


(27) 


(28) 


[k-2{hlWN)yr]Wl 
where a is a constant. On the other hand, from ([7]) it is easy to see 

t = T. (29) 

The formulas (l28]l and (f29]l constitute the desired relation between x, t and y, r. 


Examples: one-soliton solution and two-soliton solution 


The formula (1241) together with (1281) and (I29|) provides an analytic solitary solution of the CH 
equation in parametric form. As examples, we present the single-soliton solution and two-soliton 
solution in the following. 


Example 1: one-soliton solution of CH equation 

For N = 1, from (|18p we take 

fi = cosh^i, 6 = ^pi{y + cir). 


(30) 


where = Ai -|- ^, ci = f^ 2 pi_i ■ The potential function U given by the general formula (l20P 
becomes 


( 31 ) 
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From the general formula (I2ip . we obtain the one-soliton solution of ACH equation 

r = k- 2{lnfi)yr = k- ^cipjsech^ ^i. (32) 

From the general formulas ()24ll . (I28p and (|29l) . we arrive at the one-soliton solution of CH 
equation in parametric form 

u = 4[(ln/i)j^^]^ - Ak{\nfi)yr -[k- 2(ln/i)y^][ln(A: - 2{\n fi)yr)]yT 
kcip\{cip\ — 2k) sech^ 

—cipI sech^ + 2k ’ 

2/,, (^cosh^i-pisinhei)^ , 

_y I, (1 - /c^Pi)[(l +/cpi) + (1 - A:pi)e2«i] 

~ k^ "" fc3[(l- A:pi) + (1 + A;pi)e2«i] + 

t = T. 

This parametric representation of the one-soliton solution of CH equation is first found in [S] 
and latter confirmed in piinmi]. 


Example 2: two-soliton solution of CH equation 

For N = 2, we take 

1 1 

/i = cosh^i, /2 = sinh^ 2 , 6 = + cir), ^2 =-^P 2 {y + C 2 r), (34) 


where ipf = A. + jfc, e. = Ipj = X, +c, ^ The potential function U given 

by the general formula (1201) is cast into 


— -211 IT 1 _ 1 (pj - 4>f) {pI cosh^ (i+pj sinh^ ^) 

4A;2 ^ ™ 4A:2 2{p2 cosh ,^i cosh ^2 — Pi sinh sinh ^ 2 )^ 


(35) 


Based on the general formula (f2n) . the two-soliton solution of ACH equation reads 


r = k- 2 (lnW.) - h ~ pD^pI^^ ^osh^ gi + pfci sinh^ 6 ) 

2 {p2 cosh ^ 1 cosh ^2 — Pi sinh ^1 sinh ^2 ) ^ 


^ A;-r. 


(36) 
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From the general formulas (|24l) . ()28p and (j29p . we arrive at the two-soliton solution of CH 
equation in parametric form 

U = 4:[{\uW2)yr? - Ak{\nW2)yr -[k- 2{\YiW2)yr]Mk - 2{lnW2)yr)]yr 

= _2kr + ryr + j^, 
k — L 

X = In -- ^^^K + a 

[k-2{lnW2)yr]Wi 

_ [Pijpl - ^)sinh^i sinh^2 - P2{Pi - ^)cosh^i coshg2 + \{Pi -P2)coshgi sinhg2]^ 

k 16k{p2 cosh ^1 cosh ^2 ~ Pi sinh^i sinh^ 2 )^ ~ 8(^1 — Pi)(P 2'^2 cosh^ + pfci sinh^ ^ 2 ) 

t = T, 

where the function F is given in the formula ()36p . By a reformulation, this solution can be made 
to agree with the two-soliton solutions reported in 


4 Application to solving MCH equation 

In this section, we shall apply our procedure to calculate the multi-soliton solution of the MCH 
equation ([5]). Similar with the CH equation case, we first map the MCH equation onto NKdV 
equation by a reciprocal transformation. 


A Reciprocal transformation for MCH equation 


The MCH equation possesses the following Lax pair with the spatial part m 

I V’x = + ^Xrn4>, 

\ 4>x = 

and the temporal part 

I V’t = [A“^ + k{u‘^ - ul)]!!; - - Ux) + ^Xm{u‘^ - ul)]4>, 

(j)t = [A“^(tt + Ux) + ^Am(u^ - u^)]V’ — [A“^ + - n^)]</>, 

where A is a spectral parameter. 

We define a reciprocal transformation (x, t) ^ {y, r) by the relation [25l [26] 


Then we have 


dy = mdx — m{u^ — Ux)dt, dr = dt. 


d d d 

dx dy ’ dt 


A 


m{u^ 



(38) 


(39) 


(40) 


(37) 


(41) 
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Under the reciprocal transformation ()40ll . the MCH equation (l2|) is transformed to 


rur + 2m^Uy = 0, (42) 

and the potential u can be expressed in terms of m as 

u = m + l;m{—)ry. (43) 

2 m ^ 

We call the coupled system (|42]l and ()43]) as associated MCH (AMCH) equation [26]. 

By the reciprocal transformation (14011 . the spatial part of the spectral problem (|38ll is cast 
into a Schrodinger spectral problem 

4>yy = + (44) 


with the potential 


U = 


1 




4m^ 2im? 

Moreover, the temporal part of the spectral problem ()39p is cast into 

't’-r — 


(45) 


(46) 


with 


V = —2{u + Ux) = —2{u + muy). 


(47) 


The compatibility condition of (1441) and (|46p gives rise to 

{ Ut = —2^y, 

\ Vyyy-2VUy-WyU = 0, 

which is nothing but a negative order KdV equation. 

Let p = ^. From (1441) , we can easily verify that p satisfies the Riccati equation 


(48) 


Py - - — + U - p^. 


Inserting (HSl) into ([M]), we can conclude that 


m = 


2 p 


4'{y,T, A) 


A=0 


2cj)y{y,T, A) 


A=0 


(49) 


(50) 


The formula (l50P will be used to calculate the multi-soliton solution of the AMCH equation. 

Prepared with the above results, we now can execute our procedure to compute the soliton 
solutions of MCH equation. 
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Step 1: Darboux transformation of the NKdV eqnation 

We start with a seed solution V = —2k and U = ^ of NKdV equation (I48p . Similar with the 
step 1 in the above section, by applying the Darboux transformation, the iV-soliton solution of 
NKdV equation (148]) can be expressed as 




V = -2A; + 4 


dyr 


lnW„, 


where Wn is the Wronskian determinant of N functions fi, f 2 ,' 
defined by the formulas 


fj = 


(51) 


The functions /,■ are 


cosh^j, 

j = 2l- 1, 

sinh^j , 

II 

1 \2 _L 1 

li? 

> h - 


(52) 


equation 


with A = 0 and U chosen as in (j51|) are 
= 


_ A Wn 


(/i;/2r" i/iv,e 2 

W{flj2,-,fN) 


_ W(fi,f2,--- ,fN,e A Wn 


Wn' 

iW 

Wn- 


(53) 


Their asymptotic properties are 


4'i 

4>i 


^ 1 

~ 2 ^ n (^ - f )e^^ 

j-1 


4>2 


N 

n 




N 


02,2/ n(-^-f)e 


i-1 


2hy 


(54) 


Step 2: the solution m(j/, r) of MCH equation 

The formula (1501) provides a way to generate the solution m of the AMCH equation from 
(j){y,T, A)|a=o- Guided by the asymptotic behaviours (fM|) . the (j) corresponding to the boundary 
condition lim m{y,T) = k must be dicpi, where di is a constant and 4>i is given by (1531) . Thus 

lyl^cxD 

the multi-soliton solution of AMCH equation (1421) reads 


m = 


01 


WnWn 


(55) 


‘24’!,y 2 {W^Wi^f^y — WN,yW]\f) 

Furthermore, substituting (I55p into (j43p . we arrive at an analytic solution of the MCH equation 
expressed in the variables y and r 

iy3 _ 


u = 


2 {WNWpy - WpyWN) WjjWl 


(56) 
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where 

Fn =W%{2WNWN,yWN,yr " 2Wl^yWN,r + WNWN,rWN,yy - WjjWN,yyr) 

- W%{2WNWN,yWN,yr " 2Wl^yWN,r + WNWN,rWN,yy - WlWN,yyr). 


(57) 


Step 3: the relation between x, t and y, r 


In this step, we hnish the coordinate transformation between x, t and y, r. When A = 0, the 
second equation in (f38]l becomes 

4>X = (58) 

Thus a solution may be chosen as 4> = e 2 . The solution </> of equation (Hil) with A = 0 and U 
chosen as in (EU is a linear combination of the two fundamental solutions (j)i and (/>2 given by 
(153|) . The asymptotic behaviours ([Mil imply that the solution corresponding to (f) = e 2 must be 
di4>i, namely, 


X 

62 


di4>i 


di 


Wn 

Wn' 


where di is an arbitrary constant. Thus we arrive at the following relation 

1 

X = In + a. 


W: 


N 


where a is a constant. In addition, from (I40p it is easy to conclude that 


(59) 


(60) 


t = T. (61) 

The formulas (|56|) . (f60]l and (1611) provide an analytic solitary solution of the MCH equation 
in parametric form. As examples, we present the single-soliton wave and two-soliton wave in 
the following. 


Example 3: one-soliton solution of MCH equation 

For N = 1, from ([5^ we take 

/i = cosh^i, = ^pi{y + cir), (62) 

with jpI = —\\i + Cl = f, 2 pi_i ■ From the general formula (fSTI) . the one-soliton solution of 
NKdV equation (HHI) becomes 

^ = 4 ^ “ = W~ 

V = —2k -|- 4(ln fi)yT = —2/c -|- cipf sech^ ^ 1 . 


( 63 ) 
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Based on the general formula ()55p . the one-soliton solution of AMCH equation reads 


m = k + 


Ak^Pi 


(64) 


(1 — kpi)e^^^ + (1 + kpi)e~‘^^^ + 2 — Ak'^p'l 
From the general formulas (1561) . (1601) and (|61l) . we arrive at the one-soliton solution of MCH 
equation in parametric form 

Aki^pl [(1 — kpi)e^^^ + (1 + kpi)e~‘^^^ + 2(1 — 


u = k + 


(1 — kP‘p\) [(1 — kpi)e^^^ + (1 + kpi)e + 2]^ 


y {1 - kpi)e^^ + {I + kpi)e 


X = - + ln[ 


t = T. 


-a 


2fc(e?i + e-?i) 


-]^ + a, 


(65) 


Example 4: two-soliton solution of MCH equation 

For N = 2, we take 

/i=cosh^i, /2 = sinh^ 2 , 6 = + Cir), 6 = ^P 2 (y + C 2 r), 


( 66 ) 


where {pj = -jXj + ^, ci = \pI = -\\l + ^, C 2 = p^zT- It follows from the 

general formula (|5ip that the two-soliton solution of NKdV equation (1481) reads 


(67) 


Tj on,, w ^ (pj - Pi){pI cosh^ gi + pI sinh^ 6) 

AkX' ^ Ak? 2(p 2 cosh^i cosh^2 — Pi sinh^i sinh^2)^ ’ 


V =-2k + 4(ln W2)yr = -2k + 


(pI “ Pi)(P2C2 cosh^ 6 + Pi Cl sinh^ ^ 2 ) 


{p 2 cosh ^1 cosh ^2 — Pi sinh sinh ^ 2 )^ 
Based on the general formula ()55p . the two-soliton solution of AMCH equation reads 

VF 2 IH 2 


m = 


2 ( 1 H 2 IH 2 ,^ - IF 2 , 3 , 1 ^ 2 ) ’ 


( 68 ) 


where 


II 2 = 2 (P 2 cosh cosh ^2 - Pi sinh sinh ^2 ), 

W 2 = ^e^\pi{pl - ■^)sinh^i sinh^2 -P2(p? - ■^)cosh^i cosh^2 + ^(p? -P2)cosh^i sinh^2]- 

From the general formulas (I56p . (1601) and (|6ip . we arrive at the two-soliton solution of MCH 
equation in parametric form 

VFl 1^2^ - F 2 


(69) 


u = 


2 {W2W2,y - W2,yW2) W^Wi ’ 


1 1^2 
t = T, 


( 70 ) 


where the expressions kF 2 , IH 2 and F 2 are given by (f69p and (l57p . 
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5 Conclusions and discussions 

In the paper, we have proposed a simple approach to calculate multi-soliton solutions of Camassa- 
Holm equation and modified Camassa-Holm equation with aid of Darboux transformation. We 
first mapped the Camassa-Holm equation and modified Camassa-Holm equation to a negative 
order KdV equation by virtue of the reciprocal transformation. Then by applying the Dar¬ 
boux transformation to the negative order KdV equation and by inverting the reciprocal trans¬ 
formation, we arrived at the multi-soliton solutions of Camassa-Holm equation and modified 
Camassa-Holm equation. 

We believe that our approach can be applied to calculate the multi-soliton solutions of other 
Camassa-Holm type equations, for example, the Degasperis-Procesi equation m, the Novikov’s 
cubic nonlinear equation [281 129] . the two-component Camassa-Holm equations presented in 
[snj-iM]- We shall study them elsewhere. 
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